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CN I Li is studied within a three-body model Li+n + n where the core is allowed to be 

^ ' deformed and/or excite. In particular, we include reorientation couplings and couplings 
between the two bound states of ^Li. Contrary to the other examples studied within 

^ . this model, we find that core excitation does not affect the structure of ^^Li significantly. 

Tjj- ! Reorientation couplings of the deformed ^Li can change the ground state of ^^Li from a 

O I predominantly two neutron (si/2)^ configuration into a (^1/2)^- In addition, we see no 

' evidence for the existence of significant d-wave strength in its ground state, as opposed 

^ ■ to the prediction by shell model. A comparison with shell model is presented. 

1. Introduction 

^ , The prima donna of all halo nuclei is ^^Li: it offers an exciting experimental challenge 
pH \ but also questions the basic concepts of a large variety of structure theories. Many exper- 
iments have been performed over the last decade, since the pioneering total reaction cross 
section experiment at Berkeley POj. These include Coulomb dissociation experiments |2|3j . 
mass measurements |4|5j . momentum distributions of the ^Li coming from the two neutron 
removal [5I7i8j and of ^°Li coming from the one neutron removal |9llUj . proton inelastic 
scattering [TT], as well as beta decay studies |12|13|14j . We now know the system is bound 
by 295 ± 35 keV against three-body breakup |3] with a momentum distribution around 
40 MeV/c [7j (see Addendum at the end of the paper). The ground state structure con- 
tains two valence neutrons in and (pi/2)^ configurations with approximately equal 
probability |lU|12j . The r.m.s. matter radius of the system is 3.55 ± 0.10 fm 15 . 

One cannot understand the Borromean nucleus ^^Li without a good description of the 
n — core subsystem ^°Li. Experimental programs are aware of this fact but, as this nucleus 
is particle unbound, it presents a further technical obstacle. The information on ^°Li is 
summarized in ^Hl- Therein, it is possible to see the large number of experiments that 
have been performed to measure ^°Li's spectrum, but also the contradictory energy, parity 
and spin assignments made. Overall one can conclude the following: (i) the ground state 
of ^''Li contains a valence neutron outside the P3/2 subshell in a Si/2 state, at approximately 
50 keV or below; (ii) there is a state around 0.4 MeV that can be understood as a neutron 
single particle excitation into the pi/2 orbital; (iii) there is no clear evidence for a d-state 
below 3 MeV. 
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In the early days, three-body models of two neutrons and an inert ^Li core were devel- 
oped to describe properties of ^^Li .17|18|19|2(Jj . At that time nothing was known about 
^°Li, and theorists could play the game of adjusting freely the n— ^Li interaction in order 
to produce a sensible ^^Li ground state |18ll9l2lj . Models produced states where the 
two neutrons were in (2si/2)^ and in states coupled to the ground state of the 

core ^Li. The relative strength of these two components varied significantly depending on 
the n — core effective interaction. In |21j a three-body force was introduced in addition 
to the two-body n — n and n — core interactions. In |17|2Uj a density dependent n — n 
delta force was used and emphasis was given to the importance of pairing. Three-body 
inert core models have been expanded to generate three-body continuum states and 
the complexity of these three-body scattering states was analyzed within the context of 
proton inelastic scattering. 

Microscopic calculations in the early nineties were unable to reproduce a realistic bind- 
ing energy for ^^Li without artificially renormalizing the interactions, namely forcing the 
Pi/2 orbital to be closer to the continuum threshold (e.g. |23|24j ). As pairing effects had 
been identified to be crucial ^Tj, by introducing a phenomenological force in the pairing 
channel, a self consistent description of the Li isotopic chain became possible within the 
relativistic Hartree-Bogolyubov framework (e.g. j2S])- Effective interactions valid near 
the driplines have meanwhile been developed in shell model j^E]- Nevertheless, config- 
uration mixing, required to produce a realistic ground state for ^^Li, is still introduced 
by hand. It is only recently that ab-initio methods have reached mass A=ll |27|28j . 
Although the general spectra look promising, there is a clear limitation on the accuracy 
of the energy (and consequently of the two neutron separation energy) and much more 
effort is necessary to obtain reliable spectroscopic information. 

There are still some open questions regarding ^^Li, even when considering the ground 
state only. Should the quadrupole deformation of ^Li play a role in the structure of ^^Li? 
Assuming the low lying states of ^°Li are mostly neutron single particle excitations, then 
we should find a d-wave around 2 MeV excitation energy, as in the case of ^^Be. Even if 
experiment has been inconclusive in this respect, if such a state were to exist, what would 
be the implications for the ground state of ^^Li? Would one expect to also have a strong 
(^5/2)^ configuration in the ground state of ^^Li as in the ground state of ^^Be? 

On the issue of core deformation and/or excitation, ^^Li is considered within the Nilsson 
model 29] . The core ^Li has an uncoupled proton that can be either in a ii' = 3/2~ or 
K = l/2~ orbital. In that simple picture, the extra two neutrons in ^^Li, occupy a 
different K-orbit and would not be strongly affected by the orientation of the deformed 
axis. In addition, a shell model calculation using two independent frequencies (one for 
the core, the other for the halo neutrons) (201 shows that core polarization due to the halo 
neutrons is very small. 

The three-body model with core excitation j31l32l33j allows us to address these issues 
keeping all the three-body dynamics. This model was applied to the Be isotopes ^^Be 
and ^^Be. The coupled channel description of the states in ^^Be was produced assuming 
a ^"Be+n structure where the 2"*" excited state of ^°Be was explicitly included in the 
model space. A rotational model was used for the core but the strength of the coupling 
between the two ^"^Be states was taken from the measured B{E2) value. An effective 
n + core interaction for ^^Be could be made to reproduce the parity inversion of the 1/2+ 



Figure 1. Representation of the rel- 
evant coordinates in the modeL 



and 1/2^ bound states, the spectroscopy of these states, as well as the d-wave resonance. 
Based on such a realistic description of the ^°Be+n subsystem [HI] three-body calculations 
were performed using hyper-spherical coordinates and a hyper-spherical harmonic basis 
|22|- In this formulation, the Faddeev equations reduce to a set of hyper-radial coupled 
channel equations. Antisymmetrization was taken into account in an approximate way 
by projecting out the lowest ^°Be+n states with large overlap with the occupied lsi/2 
and lp3/2- Our model was able to predict most features of the ^^Be system In 
particular, the ground state was predicted to contain roughly one third of core excitation 
with neutrons in the (^5/2)^- 

If only the neutrons play a role in the low lying spectra of ^'^Li and ^^Li, one expects 
similarities between the spectra of these nuclei and those of ^^Be and ^^Be. It is then 
reasonable to think of ^°Li as ^Li+n and of ^^Li as ^Li+n+n. The main difference between 
the three-body model of ^Li+n-|-n and that of ^°Be-|-n+n is the spin of the ground state 
of the core. Exactly for this reason, computational needs for ^^Li are much larger than 
for ^^Be. Also, the first excited state of ^''Be (2+) couples the valence neutrons to new 
orbitals whereas the ground state of ^Li (3/2~) enables a larger variety of orbitals than its 
excited state l/2~. Note that, in ^Li, there are reorientation couplings just for the ground 
state, which are not present in ^°Be. Are these differences important when modelling 
^^Li? Now that computational power is very much increased, it is timely to revisit the 
problem under the constraints provided by all the new ^°Li data. 



2. Ingredients for the calculations 

The Hamiltonian in the three-body model with core excitation can be written as: 

H^B =TR + Tr + Kn(r1 + Kc(i?l, + Kc(i?2, + (l) 

where the coordinates are represented in Fig. ^and ^ is an internal degree of freedom of the 
core. These coordinates are transformed into the hyper-spherical coordinates and through 
an appropriate expansion of the wavefunction, the three-body Schrodinger equation is 
reduced to a set of coupled channel equations dependent on hyper-radius p and hyper- 
momentum K. The truncation of the basis is done through the quantum number K. 
Other details of the calculations can be found in ||32p35i| . 
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1/2+ 


0.0 


1/2- 


0.32 


5/2+(?) 


1.78 


3/2-(?) 


2.69 



Table 1 



Low lying spectrum of ^^Be [34J. Excitation energies in MeV. 



Here we discuss the physical inputs: the interactions and the model for the core. The 
interaction between the two neutrons Vnn is usually taken from a parameterization of 
the low energy nucleon-nucleon phase shifts, which are well understood. We use GPT 
interaction I^HI, neglecting the LL term. Most ambiguities reside in pinning down the 
effective interaction between the core ^Li and each neutron. We will use the notation 
nlj (g) J to denote a specific n — core configuration. For example lpi/2 ® 3/2~ corresponds 
to a valence lpi/2 neutron coupled to the 3/2~ ^Li ground state. A Woods-Saxon form 
plus spin-orbit is fitted to reproduce three main structures in ^°Li, namely the virtual state 
at -f50 keV with configuration 2si/2(8>3/2~, the resonance at +400 keV with configuration 
lpi/2 ® 3/2^ and a resonance at +3.4 MeV with configuration 1^5/2 ® 3/2^. Due to the 
lack of experimental data, the position of d-resonance in ^°Li is not known. By comparison 
with ^^Be, with the same number of neutrons, one would expect a resonance close to 2 
MeV excitation (see Table Q). We place it at +3.4 MeV: if the d-wave resonance were 
at much higher energy, the binding energy of ^^Li would not be reproduced; if it were 
at much lower energy it would become more bound than a p-wave resonance, for large 
deformation. Also we consider that the 1^3/2 orbital is at —4.1 MeV, to match the neutron 
separation energy of ^Li. 
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Figure 2. Radii of the fitted n-^Li interaction as a 
function of deformation. 
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Figure 3. Depths of the fitted n-^Li interaction as a 
function of deformation. 



When there is deformation and core excitation, l,j are no longer good quantum num- 
bers, and the nuclear state J'^ contains a superposition of different nlj / configurations. 
Nevertheless, and for simplicity, we will refer to any multi-component state J'^, by the 
component into which the state collapses in the limit of no couplings. 

The form for the n — core interaction is: 



VL f h Y ^ K„ d 



■ R.i-Rs 



Quadrupole deformation (32 serves as the collective degree of freedom for the core in- 
troduced in EqUJ The central Woods-Saxon part depends on the deformation through 
the radius R{9,(f)) = + /?2 ^2o(^;0)); where 9 and (f) are spherical angles in a rest 

frame of the core. The spin-orbit term is left undeformed. In addition, the core is treated 
within a rotational model, parameterized in terms of (32- This is an oversimplification 
of the structure of ^Li, and thus we will not impose that the strength of the couplings 
be determined by the B{E2), which in any case is not known. Volume conservation, 
however, is imposed on the central n — core interaction so when (52 7^ 0, the radius 
parameter R^g is adjusted accordingly. We take a standard value R^s = 2.60 fm when 
/92 = and its dependence for 7^ is shown in Fig. |21 Radius Rgo was made equal to 
Ryus at any deformation. The corresponding diffusenesses are fixed to the standard value 
Qws = dso = 0.65 fm. 

The depth of the central part of the interaction between the core and the neutron 
depends on their relative angular momentum /. We consider and V^^ for / = and 
1, respectively and the same V^g for all partial waves with / > 2. The depth Vgo of the 
spin-orbit term is /-independent. As mentioned before, the depths of n — core potentials 
were fitted in order to reproduce the low lying continuum structure of ^°Li as well as the 
position of bound IP3/2 neutron orbital. 
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Figure 4. Deeply bound states in ^°Li as a function of 
deformation. The legend refers to nlj valence neutron 
orbitals (coupled to the 3/2~ ground state of ^Li) and 
total spin of ^°Li. 



In principle, due to the low excitation energy in ^Li (2.69 MeV), we should include 
the 1/2" state in our core model. Such calculations were performed but results were not 
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Figure 5. Continuum states in ^"^Li, as a function of 
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orbitals (coupled to the 3/2~ ground state of ^Li) and 
total spin of "'^'^Li. Both s-wavc states refer to virtual 
states while the other states correspond to real reso- 
nances. 



7 



1.2 

0.8 

> 

S 0-4 



w 

0.0 
-0.4 

4 8 12 16 20 24 

K 

Figure 6. Convergence behaviour for ^^Li as 
a function of the hyper-momentum K. The 
hues are an exponential fit. 

sensitive to the inclusion of this additional state. As mentioned in Section ^ and due to 
its low spin value, there are no new orbitals brought into the configuration space when 
including this state. For this reason in ^^Li reorientation is significant but core excitation is 
not. Note that the model space when both core states (3/2~;l/2~) are included increases 
considerably. We will present results for the case where only the ground state 3/2" of ^Li 
is taken into account. 

In addition, the deformation parameters of ^Li can also take negative values (if the 
core is oblate). We found that the quadrupole force for oblate shapes, produces more 
repulsion when compared to the prolate case, and therefore less binding energy. Based on 
the known positive quadrupole moment of ^Li, we expect this to also be positive for ^Li. 
Therefore, we consider always positive deformation for all results discussed in this work. 

In case of an undeformed ^Li, all J'" states in ^*^Li originating from a given nlj neutron 
coupled to the core are degenerate. The degeneracy is, however, removed as soon as the 
spherical symmetry is broken by non-zero deformation. In the deformed cases, fitting 
potentials means to adjust their depths such that the centroids of lp3/2, 2si/2, lpi/2 and 
1^5/2 orbitals are kept at -4.1 MeV, +50 keV, +400 keV and +3.4 MeV, respectively. 
Again here, the 1^5/2 centroid at 3.4 MeV ensures that even when (32 = 0.7 the splitting 
does not produce inversion of the p-state and the d-state. The variation of potential 
depths with deformation is shown in Fig. |21 As in ^^Be, the spin orbit strength needed 
is large [HI]. The corresponding two-body bound states and the lowest resonances in ^"^Li 
are shown in Fig. [Hand Fig. |S[ respectively. 

Finally, in order to approximately satisfy the exclusion principle, the bound lsi/2 and 
lp3/2 neutron orbitals are projected out of the model space before diagonalization 

3. Results as a function of deformation 

Calculations were performed using EFADD ji37j. Our model space consists of 18 La- 
guerre polynomials for the hyper-radial part, reaching out to pmax = 20 fm, and 30 Jacobi 
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polynomials for the hyper-angular part. Calculations were performed up to Kmax = 22. 

In Fig. El we show the convergence of the three-body binding energy of ^^Li, for = 
0.0, 0.3 and 0.6, as a function of hyper-momentum. On one hand, the convergence exhibits 
the well known exponential dependence when Kmax ^ 12 for all deformations studied. On 
the other hand, in all cases, the convergence rate is very slow, much slower than in the 
^^Be case (see Fig. 6 of ^32 ). Moreover, the convergence rate decreases with increasing 
deformation. It is thus necessary to use the extrapolated energy value for an accurate 
binding energy. We show results up to = 0.7 (this value is already unrealistically large). 

Results for the three-body energy of ^^Li are presented in Fig. [7| We plot both the 
values for Kmax = 22 and those obtained through the extrapolation in K. Contrary to 



9 



1.0 



0.8 



0.6 - 

0.4 - 

0.2 - 

0.0 - 



0.0 



0.2 



0.4 



0.6 



Figure 9. Structure of the ground state of 
^^Li as a function of deformation. 



the case of ^^Be where the energy gain was large, here a very small additional binding is 
obtained from the quadrupole coupling, and instead, as the deformation becomes large, 
the system becomes less bound. The r.m.s. radii of the corresponding wavefunctions are 
shown in Fig. |H| Here we have used 2.32 fm for the r.m.s. radius {Rrms) of the ^Li 
core |HH1- These radii were obtained with Kmax = 22, but the Rrms variation between 
Kmax = 20 and K^ax = 22 is less than 1%. Our prediction for the charged radius obtained 
with j32 = 0.3 is Rch = 2.49 fm, in agreement with the very recent measurement 
The experimental two-neutron separation energy and the r.m.s. matter radius impose 
a constraint on /?2 ^ 0.3. Note that the extrapolation was only done for extracting 
converged energies, since it is the only observable that has a well established exponential 
dependence on Kmax- 

We now look at the structure of the state. In Fig. IHlwe show the probabilities associated 
with the three main components of the system, namely (si/2)^, and (^5/2)^. If P2 

is small, the ground state of ^^Li is almost 60% (pi/2)^; for large f32, it becomes more than 
80% (si/2)^- The region around P2 = 0.3 corresponds exactly to the transition between 
these two configurations and both are populated with equal probability. Regardless of 
the deformation, the (^5/2)^ configuration is always small (< 7%). This result is in 
contrast with the three-body calculations for ^^Be in [221 where the (6^5/2)^ configuration 
was f» 30%. 

Finally, one has to realize that there is no unique parameterization for the effective 
n — core interaction. We have convinced ourselves that the features shown here for the 
structure of this three-body system, do not result from a specific parameterization. Other 
parameterizations, using other radii or spin-orbit parameters, produce exactly the same 
features. These emerge from the constraints of continuum states in ^°Li. Variations on 
these constraints are considered at the end of the next section. 
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3-body 


shell model 


experiment 
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1.82 


2.11 
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0.37 


2.24 
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2.62 


2.73 



Table 2 



Spectrum for ^^Be: comparing three-body results with shell model and experiment |33j . 
Excitation energies in MeV. 

4. Comparison with shell model 

One may expect ^^Li to have a similar structure to ^^Be. Three-body calculations for 
^^Be have been compared with data and shell model in the p-sd shell using the WBP 
interaction. The results for the first low lying states of ^^Be from are presented in 
Table 121 As can be seen, both models are able to reproduce the first 0"*", 2+ and 1" states 
very well. Not shown, the two models also agree on the structure of the ground state of 
^^Be. A three-body versus shell-model comparison of spectra for ^^Li is not easy because 
^^Li excited states are two-particle unbound. 

For the shell-model structure of ^^Li we use the WBP Hamiltonian of [IHI- The WBP 
Hamiltonian is applicable to pure Nhv model-space truncations and is determined from a 
least squares fit of potential-model interaction strengths to experimental states which are 
relatively pure 0^, Ifiw and 2fvjj in structure. Application of this model to ^^Be and ^^Li 
gives Qhhj and 2fujj states which are degenerate within 100 keV. This implies wavefunctions 
for both i^Be and ^^Li which are 50% (s)^(p)^"^ (O^cu) plus 50% {sY{p)'^-^{sdf {2hLu) 
and other small amplitudes. 

In shell model, the Ohw ^^Li wavefunction must have 2 neutrons in the pi/2 orbit while 
for the OfUj ^Li it predicts 0.79 neutrons in the pi/2 orbit. This means that ^Li is not a 
closed core, and consequently, the Ohu ^^Li wavefunction does not look like Oliw ^Li plus 
two pi/2 neutrons. The excitation energy of the l/2~ ^Li state predicted by shell model 
is lower than the experimental value (1.395 MeV to be compared with the experimental 
value of 2.69 MeV). The 2tiiJ ^^Li wavefunction has the structure of two-neutrons in the 
sd-shell outside the ^Li core, with occupancies of 50%, 29% and 18% for si/2, d^/2 and 
(^3/2 respectively. 

The Qhjjj (negative parity) and ITujj (positive parity) spectrum for ^°Li with WBP, up 
to 4 MeV, is presented in Table El of which the first doublet 2~,1~ corresponds to Si/2 and 
the second doublet 1"*",2+ to pi/2 strength. The lowest state coming from the (i5/2-wave 
quadruplet is at i?^.(4^) = 2.05 MeV, much lower than the energy we considered in the 
three-body calculations (3.4 MeV). 

For the ground state of ^^Li, shell model predicts 15% d-wave. This 15% d-component 
in ^^Li comes mainly from the low lying d- waves present in the ^°Li spectrum. 

Also relevant for the discussion is the overlap (^^lA{g.s.)\{a^ ,a^)L\^lA (g.s)) which is 
predicted to be very strong and nearly exhausts all the spectroscopic strength. This is a 
clear indication that the first excited state of ^Li is not important in the description of 
the ground state of ^^Li. 

It is important to understand the implications of the ^°Li structure on ^^Li. Albeit the 
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Table 3 



Shell model spectrum for ^°Li. Excitation energies in MeV. 

large number of experiments, a close study of [16] raises questions about the precision with 
which states in ^°Li are known. We have explored the possibility of different assumptions 
for the states to which we fit the n — core interaction, namely Si/2 ® 3/2~ at +50 keV; 
Pi/2 ® 3/2^ at +400 keV and ^5/2 $5 3/2^ at +3.4 MeV. Of these, the least uncertain is 
the ^°Li ground state. We have checked that the main features of the present work are 
not changed by moving the p-wave to 500 keV. More important is the uncertainty in the 
location of the d- resonance. There is no clear experimental evidence for a d-wave at +3.4 
MeV, or at any lower energy. If it is a broad resonance, or it is superposed by other states 
belonging to the doublets, it could be much harder to observe experimentally. The width 
of the d-wave resonance increases rapidly with energy. If a ^°Li d-state were at 1.2 MeV 
as in ^^Be, the width would be 150 keV. The ^''Li = 2.05 MeV predicted by shell model 
has r = 420 keV and if it were shifted to E^ = 3.4 MeV as in the three-body model, the 
width would become F = 1.2 MeV. Such a broad state would be very hard to measure. 

As seen above, shell model predicts a 4~ d-wave state at around 2 MeV. What happens 
in the three-body model if the d-wave is pushed down? We refitted the n—core interaction 
with j32 = 0.3, fixing the centroid of the d-resonance at 2.5 MeV. The quadruplet states 
(l-,2-,3-,4-) move from (2.67; 3.50; 4.22; 3.21) to (1.79; 2.53; 3.32; 2.36) (MeV), 
respectively. An immediate consequence is the additional binding in the ^^Li system 
(fa 150 keV). However this extra d-wave attraction is not sufficient to change the structure 
of ^^Li's ground state, which remains essentially (si/2)^ and {pi/2)'^, with (^5/2)^ only up 
to ^ 10%. 

5. Conclusions 

We perform three-body calculations for ^^Li including deformation and excitation of 
^Li. We find that reorientation effects of the core can account for the known configuration 
admixture of s-waves and p-waves in the ground state of ^^Li. With a simple rotational 
model for the core, it is possible to obtain an effective interaction for n— ^Li, fitted to 
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the properties of ^°Li, that produces the two neutron binding energy, the r.m.s matter 
radius, the charged radius, and the structure of ^^Li all consistent with experiment. Core 
excitation in the ground state of this three-body system is found to be unimportant. This 
result is confirmed by shell model calculations. Finally the d-wave strength is predicted to 
be very small (~ 7%) in the three-body model. This is in disagreement with shell model 
predictions that show roughly 30%, 50%, 20% for (si/2)^, (^1/2)^5 and cP components. The 
main reason for this discrepancy comes from the corresponding states in the subsystem 
^°Li. So far, experiments have not been able to make a clear statement about d-states 
in ^°Li. In our n — core effective interaction we have assumed that such states would 
be above 3 MeV. However, shell model produces d-states that are at much lower energy 
(around 2 MeV), similar to ^^Be. Resolving experimentally the d-states in ^''Li will settle 
once and for all the structure of the ground state of ^^Li. This experiment should be done 
with a reaction starting from ^Li rather than knock-out from ^^Li since there is not much 
d-wave in ^^Li to start with. One possibility would be to repeat the ^Li(d,p)^°Li |1T] with 
higher beam energy. 

Addendum 

After the completion of this work we learned about the recent high precision mass 
measurement of ^^Li: 5*2™ = 376 ± 5 keV Even though the experimental accuracy is 
now well beyond the accuracy of the three body model, this additional binding may be 
an indication that the d-wave is below 3 MeV. 
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